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Abstract 

In the context of the AdS/CFT correspondence we discuss the gravity dual of a heavy- 
ion-like collision in a variant of M = 4 SYM. We provide a gravity dual picture of the entire 
process using a model where the scattering process creates initially a holographic shower in bulk 
AdS. The subsequent gravitational fall leads to a moving black hole that is gravity dual to the 
expanding and cooling heavy-ion fireball. The front of the fireball cools at the rate of 1/r, while 
the core cools as l/y/r from a cosmological-like argument. The cooling is faster than Bjorken 
cooling. The fireball freezes when the dual black hole background is replaced by a confining 
background through the Hawking-Page transition. 
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1 Introduction 



The AdS/CFT correspondence [1] has provided a framework for discussing a strongly coupled regime 
of gauge theories in terms of their gravity dual description. The equilibrium finite temperature 
problem using a black-hole background was discussed in [2] and the bulk thermodynamics was 
discussed in [3]. The transport coefficients [4] in this approach gave a result surprisingly close to 
what is measured in current heavy-ion collisions at RH1C. Also, the AdS/CFT provides a simple 
explanation for high energy jet quenching at RHIC [5] an issue of considerable experimental interest 
in the sQGP [6, 7, 8]. IN fact, following Fermi [9], who first suggested that the collision of strongly 
interacting matter will produce a thermal state, Landau [10] observed that the system would follow 
an adiabatic cooling path through transiting thermal states with entropy conservation. He further 
pointed out that the evolution should then be described by (ideal) hydrodynamics. Indeed, one of 
the key feature of the 'strongly interacting' Quark Gluon Plasma (sQGP) is precisely the observation 
of a hydrodynamical expansion in the form of radial and elliptic flow at RHIC. 

In [11], black hole formation in AdS space and its gauge theory dual was discussed using the 
setting put forward by Polchinski and Strassler [12]. More recently, in [13], the authors discussed 
a scenario leading to a long-lived or quasi-static plasma ball in the same setting. However, high 
energy collisions in QCD do not result in stopping of the through-going partons. 

The purpose of this paper is to provide a general scheme to address the complex issues of 
thcrmalization, entropy formation, cooling and freeze-out in a heavy-ion collision using the gravity 
dual description of the time dependent black hole formation. A brief summary of our results was 
presented in [14], whereby we use the AdS/CFT framework along the lines suggested in [5, 12, 11, 
13, 15, 16, 17, 18]. In short, it should be a process of black hole formation followed by a Hawking- 
Page transition which, from the boundary point of view, corresponds to thermalization, cooling and 
finally freezeout through a confinement-deconfinement phase transition respectively. Although the 
secondary scattering of partons at the boundary is a quantum mechanical process, its gravity dual 
is a classical one. 

QCD is asymptotically free and at very short distances the interaction cannot be strong. 
Thus Landau's scenario can only be applicable after some 'parton thermalization' time. In this 
respect, strongly coupled A/=4 SUSY YM theory is simpler since it is strong from fiat. Hence if we 
can prove that Landau hydrodynamics works in this theory, perhaps with some modifications and 
corrections, then we can hope to extend the arguments to more QCD-like theories with asymptotic 
freedom and chiral-deconfinement transitions. 

The rest of the paper is organized as follows: In section 2, we give a step- by-step gravity dual 
to a heavy ion collision, emphasizing the falling of closed strings created by the scattering before the 
large black hole creation. We argue that the black hole formation is an inevitable consequence of 
the falling in AdS space. We suggest that the cooling rate of the initial stage of the fireball relates 
directly to the falling rate of the black-hole in AdS. In section 4, we describe the late stage cooling 
of the fireball using the idea of brane cosmology. Our conclusions and discussions are in section 5. 
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In the appendices sketches of various ideas for future developments are drawn. In appendix A, we 
give a brief summary of some RHIC experimental facts for readers unfamiliar with these current 
experiments. In Appendix B, we comment on the recently suggested plasma ball. In Appendix C, 
we discuss collision geometries, entropy production, and give some estimates for physical quantities 
such as the saturation scale in terms of AdS black hole parameters. In Appendix D, we discuss 
temperature gradients as corrections to the homogeneous expansion described in section 4. 

2 RHIC collision and dual black hole formation 

Recently two of us have suggested [5] that real-time dynamics such as jet quenching in RHIC has a 
gravity dual description in the form of a gravitational wave falling on the black hole. The opacity 
length was found to be independent of the jet energy at strong coupling [5]. In a related but different 
picture, Nastase suggested that 5d black holes are formed through gravitational collisions of shock 
waves [18], following on the original work of t'Hooft in flat [19] and Giddings in AdS [11]. Aharony, 
Minwalla and Wiseman [13] suggested the black hole dual of static plasma balls in the pure SYM 
context. It would be interesting to describe the dynamical process by which such object if any is 
formed. 

Before proceeding further, we remark: 

1. A static black hole is assumed to be formed at once in the IR region [18, 13]. In heavy ion 
collisions, the fire ball takes a time (albeit short) to form. To describe formation and evolution 
process relevant to thermalization and cooling, the static approach is not appropriate. 

2. The assumption of a fixed temperature for the resulting black-hole [18] is unrealistic. In heavy- 
ion collisions there is no fixed temperature for the fireball, instead it undergoes an adiabatic 
path in the phase diagram. 

3. The heavy-ion collision involves fundamental (quark) probes which are embedded in the UV 
region [15, 16] as opposed to the adjoint (glueball) probes set in the IR region [13]. 

4. The bulk of the thermalization in heavy ion collision follows the thousands of elastic collisions 
each transferring energy of order N° as opposed to a full stopping of energy of order N% as 
in [13]. 

5. The arguments presented in [18, 13] apply equally to proton-proton collisions for which there 
is no evidence of hydrodynamics behavior, a hallmark of a large black-hole. 

The scattering process in the boundary occur with definite energy. How does this translate 
in the gravity dual space? 
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2.1 Where is the holographic image of initial scattering? 

In [12], Polchinski and Strassler argued that the gauge theory scattering amplitude is dominated by 
a contribution in the dual picture stemming from the height r scat ~ \fs. For dcfiniteness, we briefly 
review [12]. 

Consider the exclusive process 2 — > to particles. For the gauge theory momentum we 
associate the string theory momentum in bulk p M set by the height r through #1 



E>2 

yfdf? = — ?A (1) 

r 

The gauge theory amplitude A{p) at the boundary and the string theory amplitude in a flat space 
A s (p) are related by the postulated formulae 

/m+2 
drdQ 5 ^gA s {p) H &(r,fi). (2) 
i=l 

The string amplitude A s (p) fall off exponentially for small r, and the wave functions fall off at large 
r so that the maximum contribution occurs at finite height 



r SC at ~ R 2 P- (3) 

More explicitly, for to = 2, 



r sca t ~ R 2 V|t|ln(s/|t|)/(A-4) . (4) 

If r sca t is not smaller than the IR cut-off r min (the position of the IR brane) the image of the collision 
in bulk is located at a certain height in AdS rather than the bottom as is claimed in [13, 18]. The 
higher y/s the closer to the boundary # 2 . 

But how is this consistent with the fact that the AdS wavefunctions of the incoming particles 
are peaked in the IR not the UV as shown in the Fig. 1? 

*"Here we change the scale of the red shift factor from the convention of [12] for later convenience. In our convention, 
the gauge theory string tension a! — A 2 with the minimum height r m i n = ft 2 A defined by the minimum glueball 
mass A. 

# 2 The image, although localized at r sca t is not sharp unless the object has high conformal weight, for which the 

size of the holographic image along the r direction, Sr sca t, is estimated to be 

Sr scat 1 ,_. 
7=- (5) 

Tscat VA 

Therefore for high conformal weight in 4d, the holographic image is localized. To simplify the discussion, let's discuss 
the scattering in pure SYM without quarks. If we model the initial beam as a highly excited glueball, then the 
conformal weight A = X]j2=~i 2 can ^ e considered to be large to mimic a heavy ion-like object. This means that the 
holographic image of the incoming beam has a well defined height and the debris of the scattering should fall under 
AdS gravity. Below we will argue that the debris form a a receding black hole within a dynamical time of order ^qq d 
at the bottom or IR brane. 
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Figure 1: Wave functions of a few low energy glueball spectra. 



The answer lies in the measure factor. Since the discussions so far did not introduce fcrmions 
(fundamental fields), we model the heavy ion collision as a collision of two highly excited glueball 
states. For an incoming glueball with definite energy, the scalar wave function is factorized as 
$(x M ,z) = e <f){z), and the scalar field equation in 5 dimension, (D 5 — m 2 )$(a; A1 ,z) = 0, reduces 
to 

[z 5 d z z~ 3 d z - k 2 z 2 - [mRf] 4>(z) = 0, (6) 

with k 2 = k 2 — uj 2 . Therefore for a process with definite energy- momentum, the wave function 
depends only on the mass not on the energy. For a confining theory, we need to cut off the IR part 
by hand (restrict to z > z m ) or by a pertinent metric structure, hence we are interested in a wave 
function that is regular near the boundary (z <~ 0), where, the wave function behave <j) ~ z 2±I/ with 
v = (mR) 2 and i? 4 = 4irg s Na' 2 . The non-normalizable part (z 2 K u (kz), z 2 N u (kz)) should be 

interpreted as the two point function with a source at the boundary [25], so that its coefficient is the 
strength of the source. For the initial beam, it is on-shell and k 2 < is the 4 dimensional mass. The 
explicit normalizable wave function is [24] 

<j>(z) = z 2 J u (kz) for k 2 < 0. (7) 
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Figure 2: Probability distributions. Notice that z = is the boundary. For higher excitations it is more 
likely to be at UV region. 

On the other hand and between the collisions, the particles are off mass-shell (k 2 > 0)and k should 
be interpreted as the momentum transfer. The wave function is 



As the 5d mass (equivalently ~ v) increases, no qualitative change in the wave functions is 
observed except that it is slightly pushed to the IR region(larger z). On the other hand, increasing 
the 4d mass (fc 2 ) includes more nodes in the allowed region and effectively pushes the wave function 
into UV region. This "push-to-UV" effect is more dramatic if we consider the 'radial' probability 
density P(z) = ^/g\(j>(z)\ 2 . Due to the measure (y/g), the dominant peak is near the boundary rather 
than horizon. We can estimate the location of the dominant peak in terms of x„i, the first zero of 
the J„ (x) . We suggest that the location of the holographic image of the incoming glueball with mass 
M4 (= \J — k 2 := k) is given by 



which is consistent with eq.(3). 

In summary, if we model a heavy ion as a glueball with large 4d mass, the holographic 
image of the initial beam is at the height that is proportional to the mass. This is consistent with 



(j>{z) = z 2 I u (kz) for k 2 > 0. 



(8) 




(9) 
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Polchinski-Strassler's argument as detailed above. The scattering takes place at a given height since 
their initial states are localized there. 

2.2 Expansion and Thermalization 

Since RHIC data shows that the fireball after collision is a thermalized, an AdS black hole should be 
formed. But then how?. For this, we notice that the AdS gravity has anti-tidal force so that it has 
a focusing property. Namely, two vertically separated particles in AdS bulk will become closer as 
they fall. We have a gas of falling debris after collision. If one consider a local rest frame of the fluid 
(lagrangian coordinate in fluid mechanics language), the common proper time for all the particles 
can be used as a time coordinate. This is similar to the treatment of Bjorkcn in 3+1 spacetime by 
assuming LRF. As we will see below, all the falling particles released from different heights arrive 
at the bottom afther the same proper time, r = irR/2. The unavoidable consequence of this is that 
the motion of AdS fluid is like a 'cosmological' contraction leading to the final singularity and the 
dual of the fireball forms an AdS black hole. Of course, The dual of this contraction is the fireball 
expansion in the boundary. 

To be more explicit, consider a radial in-fall in AdS space: 

* 2 = (s) V -(?)* 3 - < 10 > 

For massless particle, the motion is described by a null geodesic with solution 

r = R 2 /t, (11) 

and the falling should start on the AdS boundary at t = 0, which is consistent with the picture that 
the free falling of massless particles in AdS is dual to the free expansion in the boundary whose front 
surface is expanding with light velocity [5, 20]. For massive particle it leads to 




where e = (r/R) 2 ^ denotes the energy per unit mass. The resulting motion is harmonic in proper 
time, 

r = Recos(r/R), t = R/e ■ tan(r/i?). (13) 

The period is 2ttR which is independent of the initial conditions. In case there is an IR brane, 
the initial difference in height Sr(r = 0) = RSe will be reduced to to Sr(r) — RS€COs(t/R) at the 
bottom. In terms of the boundary time t, 
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Figure 3: Holographic correspondence of the expansion in 4d and the falling in 5d. From the boundary 
point of view, the front part 'F is freely streaming while the inner part '3' sees medium effects. From the 
bulk point of view: the lower part 'F falls freely while the upper part '3' sees the AdS black hole geometry. 
Birkhoff's theorem tells that whether the inner part is really black hole or not is not an issue. Fhus the 
inner part '3' feels that it is in thermal equilibrium. 



' 7i> = R 2 /t-{R/tf/2e 2 + 0{r 5 ) 1 (14) 



V^/R) 2 + 1 

so that the initial condition dependence ( that is the e dependence) disappears rapidly as time goes 
on. We believe that this focusing effect plays an important role in the initial formation of the black 
hole geometry. So eq. (14) can be thought to describe the front surface of the fireball which is not 
equilibrated. 

After reaching bottom (IR region) the droplet will spread and flatten to make a pancake. For 
late time falling objects, such a stack of mass on the IR brane generates a black hole geometry due 
to Birkhoff's theorem. See Fig. 3 and its capture. The particles inside the front surface, experience 
the interaction of a medium and the expansion in the center of the fireball is dual to the falling of 
a particle in the AdS black hole background. 

In the next subsection, we will consider the case with fundamental fields (quarks). 

2.3 With quarks: creation of closed string 

Now if we have particles in the fundamental color representations in addition to the ones in the 
adjoint color representation, we need to introduce probe branes in bulk [21]. A heavy meson is a 
quark and antiquark connected by a string deep in AdS. The scattering of such mesons could be 
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Figure 4: At each interaction vertex of two scattering mesons a closed string must pop up. This is a unique 
feature of AdS space that does not take place easily in flat space. 



realized by moving these AdS strings, which are highly non-local objects in AdS bulk. The collision 
of such strings in the bulk may form a highly distributed object in the bulk and may not be a black 
hole initially. However, the contraction and AdS fall of these objects will give a black hole. 
The mesons in N=4 theory were studied in [17]. They are deeply bound with mass 



M~ m « . (15) 

So if we model a heavy ion by this meson, the quark mass should be taken as large. In this case, 
the holographic image of the fireball is created at a significant height (~ M q ), and falls to form a 
black hole at the bottom. In this picture we can argue that for each vertex, a closed string can be 
created to leave behind a flavor branc. 

The many parton collisions at the boundary trigger i. elastic collisions which are dual to 
massive closed strings; ii. inelastic collisions which are dual surface flips. An example of the former 
process is shown in Fig. 4. Since the minimal string is not a straight line connecting two sources on 
the boundary (infinite warping factor), the string must stretch inside the AdS space [26, 27]. As two 
mesonic composites come together, the recombination from AB + CD to AD + BC should happen 
just before B and C touch each other, since that is energetically favored. For example, when the 
separation (in boundary) of AB and BC are both L and that of BC is e, then for small enough e 
the difference of total lengths of the string is 

Iab + Icd-Iad-Ibc = -^y + or°, + ~ > °> ( 16 ) 

where c is just a constant. Thus in a hadron-hadron scattering process, recombination of the string 
must arise at the vertex (where B and C coincide) generating a closed string. This is a remarkable 
feature of AdS space with no analogue in flat geometry. Although the above example is for pure 
AdS, we expect the mechanism to be universal regardless of the geometry in the IR region if the UV 
region remains AdS. 

Each of these liberated closed strings fall in the AdS space under AdS gravity. Some of the 
the closed string states could be in a black hole state. They merge as they fall to form a larger black 
hole in the IR region as we suggested earlier. In fact, the initial colliding objects contain multitudes 
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Figure 5: Multiple interaction vertices create a shower of massive closed strings in AdS space. Some of 
them are mini-black holes. The strings flake and fall towards the AdS center like a rain-fall to form a large 
black hole at bottom. 
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Figure 6: Schematic view of wall- wall collision. The vertical coordinate is the 5th dimension r, the horizontal 
one is the longitudinal one. (a) before the collision, (b) shortly after the collision: color rearrangements and 
massive production of closed strings. They fall into AdS space; (c) A larger black hole is formed at the 
bottom. It will flatten to a pancake shape, lowering its horizon. 

of particles, and they can not be in thermal equilibrium at once. Therefore in the gravity dual what 
forms immediately after the collision is not a single big black hole but a multitudes of mini black 
holes together with other closed string states. 

The efficient creation of the particles provide a mechanism to convert the deposited collision 
kinetic energy to mass resulting in lowering the temperature scale and most of the initial energy 
is deposited as mass. In real QCD, this procedure increases the strength of the interaction by the 
running coupling. Here in N=4 SYM, there are no such effect since coupling does not run. From 
the boundary language, the increase of number density of particle and the increase of interaction 
strength is the key point to get the efficiency in the thermalization. 

The gravity dual of the collision processes is shown schematically in Fig 6. Even though 
matter is partially stopped on the boundary (UV region) , the initial entropy build up will cause the 
formation of a fireball and its expansion. Its gravity dual is a detached set of closed string states 
that fall into the AdS space. 

2.4 Final Stage: Hadronization 

As time goes on, the black-brane like object becomes thinner and eventually unstable to density 
fluctuations that lead to an instability [22]. We expect that the extremely thin black pancake will 
fragment into small pieces each of which evaporates (quantum mechanically) by Hawking radiation. 
This may be identified as the confinement phase transition. Since the resulting metric is nothing but 
the AdS metric with IR cut-off, which is a confining metric, this can be considered as the details of 
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gauge d=4 theory 



string/gravity in d=10 



Total cm energy per unit mass 

initial gluons, CGC 

Thcrmalization — > T 

the entropy 

rescattering (q — q) 

rescattering (3 — 5) 

fireball expansion before equil. 

fireball expansion after equil. 

further equilibration 

ideal hydrodynamics 

hydro with viscosity 

kinetic freezeout 

deconfincmcnt 



height of scattering brane in AdS 
Aichelburg-Sexl-type shock waves 
black hole formation with T bh = T 
area of the horizon 

production of closed string (gravitons) 
interactions between closed strings 
falling b.h. in AdS bulk 
spreading of b.h at bottom 
merger of gravitons to black hole, 
stationary black hole 
growing black hole 
cutoff of gravitons 

fragmentation of thin black-brane followed 
by hawking evaporation 



Table 1: A vocabulary of dual phenomena in gauge and gravity formulations. 

the Hawking-Page transition [2, 17]. # 3 In our picture, we do not expect a quasi static black holes 
as in [13]. This is also so in heavy-ion collisions, since not all partons are stopped at once, the initial 
expansion is one dimensional instead of three dimensional as Bjorken suggested [23]. After local 
equilibration is achieved in a heavy ion collision, the matter expands and the temperature depends 
on both the location and time. 

With these considerations, the gravity dual of the RHIC collision can be set up by considering 
the physical process together with the general dictionary of ADS/CFT listed below. 



3 The cooling of the core of the RHIC fireball 

In this section we discuss cooling and expansion of the fireball in the late stage, where the black hole 

becomes a black-brane like object which is expanding in spatial direction and lowering its horizon 

continuously under the influence of AdS gravity. The temperature decrease is interpreted as the 

increase of the distance between the probe brane and the black hole horizon. 

Since the AdS space is homogeneous, we can change the frame such that the black brane 

is fixed while the probe brane is moving in the background of the static black brane. Then the 

probe brane moves to the UV direction and therefore sees bigger scale factor of the bulk metric 

as it moves, resulting in the cosmic expansion on the brane world. This is nothing but the brane 
#a In 

pure N=4 SUSY the interaction is the same at all scales. Thus an expanding fireball of "CFT plasma" will 
never freezeout, and will expand hydrodynamically forever till zero temperature is reached. Freezeout can be reached 
in the gravity dual by switching to confining D-brane metrics, which is known as a Hawking-Pagc transition. 
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world cosmology addressed in [29, 30, 31, 32, 33, 34]. In this way, we identify the gravity dual of 
the expansion/cooling of the fireball as the cosmic gravitational expansion in the AdS-black-hole 
background. In other words, we approximate the Little Bang as the Big Bang on the probe brane. 
In the Little Bang, the temperature has space and time dependence, while in the Big Bang there 
is no spatial dependence, only time dependence. Therefore the approximation is good only for the 
center of the fireball, which is the subject of this section. 

Although the real expansion is mostly 1 dimensional, we believe that the thermally equili- 
brated expansion is 3 dimensional in nature. This is because the 1 dimensional expansion is driven 
by the ultra relativistic motion of the initial particles whose speed can not be caught up by the 
interactions. 

3.1 Big Bang on a moving brane 

We consider a class of metric given by the near horizon limit of non-extremal D p branes: 

ds 2 = g Qa dt 2 + g{r)dx p 2 + g rr (r)dr 2 + g s dfl s _ p , (17) 

where g = (r / 'R^ 7 ^/ 2 , \ 9oo \ = {r / R)^-v)' 2 {l - (b/r) 7 -P) = g- 1 and g s = ^{R/rf 7 ^/ 2 . The 
dilaton is given by 

** - (7) • (18) 

If we neglect the brane bending effect and consider the configuration of zero angular momen- 
tum of the brane around the sphere, the DBI action for the D p brane 

S P = -T P J (P+ 1 Ze-+y/-det7 a fi-T p J C p+U (19) 

can be written as 

S p = -T p J dP +1 de-*gPl 2 ^\g m \-g rr P. (20) 
Since there is no explicit time dependence 

E=p-q-L= 9 -C, (21) 

V l.9oo I - g rr r z 

with C = (r/R) 7 ~ p , is a constant of motion. Using the equation of motion 

grrr 2 + goo + g p \goo\e- 2 */(C + E) 2 = 0, (22) 
the induced metric can be written as 

2 p -20 

ds 2 = - g -0^ w dt 2 +gd X 2 . (23) 
Defining the proper (cosmic) time r by 

dr= \g o\g p/2 e-*/(C + E)dt, (24) 
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the induced metric can be written as a zero curvature Friedman-Robertson form 

ds 2 = -dr 2 + a 2 (r)dx 2 , (25) 



where a 2 = g(r(r)). The equation of motion can be rewritten in terms of a and r 



{C + E) 2 e 2 ^ 1 \ 



with g' = dg/dr. Then, the equation of motion in terms of a and r is given by 

^ = ^ 7 ~P V a 2(3-p)/(7-p) 



r E / /,• P 1 



(26) 



(27) 



i 4 ; V i? 7 ~ p a 4 , 
where we have used the fact C = (r/R) 7 ~ p = a 4 . Notice that the effect of the RR-fiux field C is to 
provide a strong enough repulsion force to cancel the confining AdS gravity. 
As a — > oo (late evolution), we have 

a(r) ^ T V-p)m-p). (28) 

The scale factor evolution a(r) captures the cooling of the fire-ball at the boundary through its 
holographic dual: 

T(a) = « (29) 

Vl5oo| °( T ) 



with 

as the black hole temperature. The local temperature is the black hole temperature observed by the 
observer in the probe brane. This is the actual temperature of the fireball. As the brane moves away 
from the black hole, the brane world (the fireball) expands and cools according to T(a) — T^ja. 
There are two interesting cases: p — 3 and p = 4. For p = 3, 

a(r)~V^, T ~ -^=. (31) 

The reason for considering p = 4 is that one of its direction (say £4) in a confining theory is 
compactificd. After the compactification the p — 4 and p = 3 are identical. Without compactification 
a w t 3 / 7 which is a stronger warping. 

This result is to be compared with cooling in D-space. Indeed, the entropy for a (perfect) gas 
is just S w T D Vd- For a relativistic d-space hydrodynamical expansion we expect Vd ~ Vo-d T d . 
For fixed entropy, the temperature falls like T as 1/T d / D . Bjorken 1-space expansion (31) corresponds 
to D = 3 and d = 1, therefore T a 1/r 1 / 3 . Fully 3-space expansion corresponds to T a 1/r. The 
AdS case with T a 1/\/t is faster than Bjorken in lspace but slower than perfect hydrodynamical 
expansion in 3space. It is like fractal with d = D/2 = 3/2. One may summarize these result by 
saying that strong interactions slow down the expansion of the fireball just as gravity does in the 
dual picture. We note that since the viscosity is quantum with rj/(S/V 3 ) a h/iir its effects are not 
present in our estimates. Their consideration follow from perturbation theory and are easily seen to 
delay the cooling. 
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3.2 confinement phase transition 

When T cools enough such that T < Aqcd, there must be a Hawking-Page transition [28] and the 
background metric is replaced by 

ds 2 = {r/Rf' 2 {-dt 2 + dx 2 + fa{r)dx\) + {R/rf' 2 {dr 2 / fa + r 2 dQ 2 4 ), (32) 

where fa = 1 — {tkkIt) 3 refers to the compactified direction. Witten [2, 17] suggested that the 
transition to this metric maybe interpreted as the confinement /deconfinement phase transition. The 
equation of motion in the new background can be calculated. Though minor, there are a few 
differences in detail of the calculation, but rather surprisingly, the final outcome is precisely the 
same with the substitution b — > Tkk- For p = 4, 



= 16° ' 



E A 2 (. r\ K 1 



+ 1 



a 4 V R 3 a 4 



(33) 



As we discussed before, The front factor 9/(16a 2 ' 3 ) disappears if X4 is compactified (which is effec- 
tively p — 3). The phase transition point in terms of the brane position occurs when the warping 
becomes dp at 

T{a F ) w T KK , (34) 
where the Kaluza-Klein temperature is given by Tkk = Sr^^/^irR 3 / 2 ). Thus 



a F = ^ = J— (35) 
Tf V r KK 

One may interpret the phase transition as hiding of the black hole horizon behind the KK singularity 
r = r^K ~ 1/Aqct>- After this phase transition, hadron creation begins which should be a dual to 
the Hawking evaporation of black-brane after Gregory-Lafflame transition as we discussed at section 
3. The fireball ultimatly freezes out when the pions decouple. 

So far the expansion is homogeneous. This is consequence of the assumption that the back- 
ground is black brane. In reality, temperature has some gradient. To discuss the spatial dependence 
as well as the time of the temperature we need the brane motion in the background produced by an 
object with finite extent like a black hole. We give a touch to this important but difficult subject in 
appendix D. 



4 Discussion 

We close by summarizing. Recent heavy ion collisions at RHIC have suggested that the released 
partonic matter produced is strongly interacting in the form of an sQGP. Two of us have argued 
recently that a good starting point for adressing key issues of the sQGP is A/*=4 SYM at strong 
coupling. In this paper we have suggested that the entirety of the RHIC collision process from the 
prompt entropy release to the freezeout stage can be mapped by duality to black hole formation and 
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evolution in AdS space. In other words RHIC little bang and cooling is dual to a cosmological big 
bang with a flying black hole as a proceed. 

We have provided a simple physical picture of black hole formation and thermalization from 
string theory point of view. We have suggested that due to the strong interaction of the fireball 
liquid, the expansion is slower than expected from the ideal gas model. The cooling of the fireball is 
\j \pr which is slower than Bjorken 3d cooling 1/r. The strong nature of the interaction slows down 
the expansion rate hence the cooling is slower than expected from the Bjorken solution. Cooling 
freezes when the background is replaced by the confining background through the Hawking-pagc 
transition. 

There is a clear distinction between coherent parton-parton scattering and incoherent macro- 
scopic (heavy-ion) collision of large number of partons. In the former, the scattering happen in IR 
region and information is conserved, which is a hallmark of quantum mechanics, while in the latter 
the scattering happens at the UV region and the information is lost #4 and entropy is generated, 
perhaps to its maximal value, typical for local thermal equilibrium. While entropy generation maybe 
traced back to the incoherence due to the many binary scattering in a RHIC heavy-ion collision, it 
is readily understood in the gravity dual description: Since the particles evaporate from the crowds, 
the (contracting) core is losing information and become a black hole. Although the lost information 
will be back to the black hole (due to the AdS gravity) and the black hole grows, the entire infor- 
mation is hidden inside the horizon. This we believe is one of the simplest explanation for entropy 
production at RHIC. 

Non-cosmological like expansions with realistic fire ball geometries on the boundary are more 
involved to analyze. We have suggested that their asymptotic stages can be mapped on black hole 
perturbation theory resulting in non-ideal hydrodynamics from conventional Einstein gravity. We 
will report on these issues and others in future. 

Note added: Since submitting our paper, a new paper by Janik and Peschanski [47] appeared. 
The cooling of the fireball is discussed using the asymptotic solution for the metric induced by the 
energy momentum stress tensor set on the probe brane. In the gravity dual, the black hole moving 
horizon reproduces Bjorkcn's scaling for a one dimensional expansion, in support to our arguments. 



* 4 A hypothetical full experiment with measuring phases of all thousands of secondaries can still in principle recover 
it. 
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Appendices 



A Elements of RHIC physics 

This discussion is intended to be elementary to shorten up the vocabulary gaps between the string 
community and the heavy ion community interested in the gauge-gravity problems through the 
AdS/CFT correspondence. 

Collision: Experimentally we use the heaviest (and fully ionized) nuclei (mostly Au 197 at 
RHIC) with as large energy per nucleon as possible (the relativistic gamma factors 7 ~ 100 in center 
of mass, to be increase further at LHC soon.) 

One may ignore the complexities of nuclear physics and QCD evolution, and focus solely on 
the partonic wave function of hadrons or nuclei before the collision. More precisely, as coherence 
is lost anyway, one needs to know the mean squared amplitudes of the pertinent harmonics of the 
comoving gluon field with the so called saturation scale Q s or equivalently the transverse density 
of partons Q 2 . At RHIC Q s is about 1.5 GeV for a typical Feynman x = 10~ 2 . It will be higher 
at LHC say Q s — 6 — 8 GeV at lower x. A model currently used to describe the low-x part of the 
nuclear wavefunction prior to the RHIC collision is the color glass condensate (CGC). It is rooted 
on a weak coupling argument in QCD contrary to what is stated in [18]. 

Equilibration: This is a transition from the CGC to thermal quarks and gluons. Solutions 
of classical Yang-Mills, both for random fields [35] and sphalerons [36] have actually produced 
thermal-looking spectra but more is to be understood, perhaps along the discussion of plasma 
instabilities [37]. 

Hydrodynamics: This is a key aspect of RHIC physics. Maintaining collective flow for 
systems containig just ~ 100 — 1000 particles is a nontrivial issue [38], and would not happen for 
usual liquids like water. Thus the matter produced at RHIC is now refered to as a strongly coupled 
quark-gluon plasma (sQGP) or liquid. Indeed it exhibits both bulk thermodynamical parameters 
and transport coefficients (viscosity) that are surprisingly close to what the AdS/CFT correpondcncc 
predicted for strongly coupled N=4 SUSY YM theory. The short time behavior of the hydrodynam- 
ical expansion is close to the 1-d Bjorken regime whereby the temperature depends only on the 
proper time r = \/t 2 — z 2 . For central collisions the expansion becomes axially symmetric before 
turning to a full 3d spherical expansion. For non-central collisions there is azimuthal anisotropy 
which is successfully described by hydrodynamics. 

Freezeouts: This corresponds to chemical and thermal freezeouts whereby the change in the 
composition is turned off (chemical) and the particles decouple (thermal) with free streaming. Both 
freezeouts follow from the same condition ^expansion = ^reaction, where we have used the covariant 
definition of the expansion rate ^expansion = d^u 1 *. 

In cosmology, the expansion is so slow that not only strong (pp) scattering survives, but even 
weak equilibrium through p + e <-> v + n docs, untill T 1 MeV. Photons freezeout at much lower 
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temperatures T 0.1 eV. At RHIC chemical freezeout corresponds to the end of particle changing 
reactions such as 2tt — > 47r, while kinetic freezeout corresponds to the last elastic collision such as 
2-7T — > 2n. Experimentally both freezeouts are reasonably well measured, the former from matter 
composition while the latter from particle spectra [39]. While the critical temperature in QCD 
T c h ~ 176 MeV is independent on the collision centrality the freezeout temperatures depend on the 
system size. For instance, the kinetic freezeout temperature T^ n does depend on the system size, 
and goes down for the largest fireballs (central collisions) to about 90 MeV. Thus the whole range 
of temperatures at RHIC is about 4-fold, from the initial Tj w 350 MeV to the kinetic freezeout 
Tkin w 90 MeV. The energy density changes by about 2 orders of magnitude. 

The main reason for the rapid freezeout of a hadronic gas is the Goldstone nature of the pions. 
The self- interaction through derivatives makes it difficult to generate soft pions. At low tempera- 
ture, the pion gas collision rates can be calculated from the leading chiral interaction (Weinberg- 
Tomozawa). Specifically, the elastic rate is [40] 



(36) 



The strong T dependence follows from dimensional arguments. The inelastic rates can be found 
in [41]. 

At RHIC detailed numerical calculations show that the proper time spent in the sQGP phase 
(T > T c ) the "mixed phase" (T rts T c ) and the hadronic phase (T < T c ) are all comparable. However 
at LHC the sQGP should dominate. For simplicity, we may ignore the complications inherent to 
the running coupling in QCD, the confinement-deconfinemcnt transition and the pion dynamics by 
restricting the discussion to the early phase of the collision dominated by the sQGP. If the latter 
phase is close to strongly coupled N=4 SUSY matter at finite temperature, as two of us discussed 
recently [5] , it is then useful to use the duality insights to bear on the bulk and kinetic properties of 
the sQGP. 



B Comments on plasma bubble 

The formation of a metastable "plasma bubble" was recently discussed by Giddings [11] and its 
boundary and slow evaporation of glueballs was discussed in [13]. In QCD the idea that near the 
phase transition there can be a near-stable fireball was discussed since 1970's, see in particular [45]. 
Although at finite N the pressure is always finite, the ratio to the energy density p/e has a minimum 
called "the softest point" [46]. If the system is produced at such conditions the produced fireball 
would be especially long-lived: see detailed hydrodynamical studies in [46]. Experimentally there 
are indirect hints that fireball lifetime is indeed maximal around collision energy y/s w 6 GeV * A, 
which unfortunately was not studied in detail yet: there are proposals to run RHIC at such a low 
energy to verify that. The gravity dual of such slow evaporation is Hawking radiation in the (usual 
not extra) spatial direction. This phenomenon is known in QCD as a "long-lived fireball" . It was 
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first related with the MIT-type bag model, in which a meta-stable zero pressure state is possible. 
It was suggested as early as 1979 as a possible explanation for why the secondaries produced in 
pp collisions at ISR have thermal looking spectra for transverse momenta without hydro-expansion 
[44]. Hydrodynamical transverse expansion was observed in heavy ion collisions in 1990's in AGS 
experiments in Brookhaven: its small magnitude was attributed to the so called "softest point" of 
the Equation of State, the minimum of p(e)/e right after the phase transition, which is close to the 
initial conditions in those collisions* 5 . 

Aharony et al [13] argued that since the critical pressure p c w AT? while the energy density 
on the plasma side of the transition is e w N%, their ratio is 1/N%, vanishing in the limit of a large 
number of colors. Furthermore, they argued that a small but nonzero critical pressure can be further 
compensated by a surface tension <r, making stable plasma drops provided that the drops are small 
enough* 6 , with a size R < Rc w a/p c . For large N c this can again be large since a ~ N%. 

One important feature of the arguments presented in [13] is that the evaporation of glueballs 
from the surface of such plasma droplets leads to a reduction of the radius R in the ratio a/p c , 
so that inside the droplet the temperature should actually increase. This is indeed qualitatively 
similar to what happens with an evaporating black hole, emitting Hawking radiation. 

While interesting, this however is a doubly tuned situation, which is rather different from 
what we expect at RHIC. In this case, there is a violent expansion and cooling of the fireball as we 
have attempted to describe in this paper. In our setting QCD matter (and its collisions) are placed 
as a test brane removed from the IR limit of AdS. The expansion and cooling are both due to a 
warping of the metric induced by a departing black hole albeit far towards the IR. 

C A few estimates 
C.l Collision geometries 

To simplify the initial geometry of the problem, imagine that colliding bodies may have infinite ex- 
tensions in some directions, with the solution naturally independent on the corresponding variables. 
Let us call the number of "non-contracted" variables d. 

The simplest geometry (i) would be a spherical collapse: One may imagine a spherical shell 
of matter collapsing into itself with an initial radial velocity v and Lorentz factor 7 = 1/ vl — v 2 . 
A fireball which is produced in this case is expanding in a spherically symmetric way, producing a 
"Little Bang" like at RHIC, only in a much simpler spherical geometry* 7 The next geometry (ii) 
to consider is a collapse of a cylindrical shell, leaving one "non-contracted" variable, d = 1. The 

# 5 Now there is a proposal under evaluation to run RHIC at very low collision energy to possibly produce such a fire- 
ball close to its critical state and study the "softest" regime of a strongly coupled quark-gluon plasma experimentally. 

# 6 This may explain why small systems, such as those produced in pp, do not show hydro-expansion while large ones 
produced in central PbPb collisions at similar energies do. 

# 7 This has been considered by one of us many years ago [43] for e+e- collisions and prior to QCD. However, due to 
asymptotic freedom this condition cannot be created experimentally: e+e- collisions in fact result in 2 jets, propagating 
from the collision points in random directions rather than a spherical expansion. 
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gravity dual to it should have a black hole with one less dimension. The geometry (iii) with d = 2, 
is a collision of two infinite 2d walls. This is close to what happens at RHIC, where the colliding 
Au nuclei are Lorentz contracted by a factor hundred into two thin pancakes. A variant of this are 
light-like wall- wall collisions that pass through each other causing surface/string rearrangements in 
the minimal impact parameter region much like the parton-parton scattering approach originally 
suggested in [15, 16]. 

The realistic # 8 case (iv) corresponding to RHIC is a collision of finite-size objects (although 
as large as practically possible). Due to relativistic boosts the nuclei get flatten in the collision 
direction a; 3 . Furthermore, for non-central collisions the overlap region is not axially symmetric but 
has an almond-like shape. Its gravity dual presumably would create a black hole with a horizon of 
some ellipsoidal shape, with different dimensions in all directions. 

C.2 Fermi-Landau model and the entropy formation 

In QCD and other asymptotically free theories we know that at small distances (close to the origin) 
the interaction is weak. In the collision the constituent partons would literally fly through each 
other. Thus the issue of entropy formation at RHIC is complex and, as one may have suspected, 
not unanimous. 

In contrast, in strongly coupled N=4 SUSY YM theory, there is no relation between the 
coupling and the scale. At strong coupling, one may think that the colliding matter is stopped, 
and that most of the entropy is produced promptly at this stage. Thus, we use for this case the 
Fermi-Landau (FL) model [9, 10] as a benchmark for further comparison. 

The main assumption of FL is that matter can be stopped in a Lorentz-contracted size 
R = Ro/"f, where i?o,7 are the original size of the colliding objects and their Lorentz factor. The 
volume in which it is supposed to happen is 

V ~ Rl/j 3 ^ (37) 

where d is the number of "non-contracted" coordinates introduced in the preceeding subsection. The 
first step of the argument is to evaluate the temperature at this stopped stage. The energy density 
is 

e = E/V ~ 7 4 ~ d ~ ~ T 4 (38) 

where the last equality is from the EoS of matter. Therefore the temperature grows with the collision 
energy as 

T ~ j 1 -^ 4 (39) 
The next step gives the amount of entropy produced: 

S « T 3 V w 7 d/4 (40) 

# 8 The gauge theory under consideration is still not QCD but a strongly coupled A/"=4 SUSY YM 
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One can see, that in the spherical collapse (i), there is no entropy growth because d = 0. The lesson 
from it is that only the cases with less trivial geometry provide some interesting predictions. 

Despite the differences between the FL model and QCD, the entropy prediction for the wall- 
on- wall case (hi), S ~ 7 1 / 2 ~ s 1 / 4 , agrees with the observed multiplicity growth quite well. We will 
return to the discussion of this point later. 

C.3 black hole formation 

For simplicity, we assume that the black-hole has been formed but that the flaking of closed strings 
is still taking place, and ask: under what conditions the flaking strings can be captured by the 
black-hole? what is the typical accumulated entropy? what is the typical time for this entropy 
formation? Most of the arguments in this subsection are heuristic. 
^From the AdS black hole metric we have 

ds 2 = -fdt 2 + f- 1 dr 2 +r 2 dfl 2 5 , with / = 1 - G 5 M/r 2 + (r/R) 2 . (41) 

where we have ignored the distorsions caused by the boundary brane on the black-hole. The horizon 
size of the black hole is 

Hence r^h = (B^G^M) 1 / 4 := b for a large black hole, while r^h = \fG^M for a small black hole. 
# 9 The temperature of the large black hole is given by Tbh = b/irR 2 , while that of the small black 
hole is Tbh = R/2irb 2 <~ I/VG5M. The large AdS black hole does not evaporate while small black 
holes can. However, the hawking temperature goes up as it evaporate while the fireball cools as it 
evolve. Therefore small black hole seems to be improper to describe the RHIC fireball. Therefore 
throughout this paper we identify the RHIC fireball with large AdS black hole. 
We can express the mass and entropy in terms of Hawking's temperature 

M = R e T*/G 5 , ! = ^^«T3. (43) 

On the other hand, using G^ 1 = G^R 5 = M®R 5 , M p = l/l s gl^ 4 , R 4 — g s N c l 4 , we can express 
Hawking's temperature in term of mass 

b 1 fM\ 1/4 

Since M w N 2 (see below) then T w N®. The time the entropy is reached corresponds to typically 
the falling time in AdS space 

r = \R (45) 

# 9 Large black hole means G5M/R 2 > 1 with / ~ (r/R) 2 (1 - b*/r 4 ), and small black hole means G5M/R 2 < 1 
with / ~ 1 - G 5 M/r 2 . 
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Which is of order N®. This is about the time it takes the final black hole to reach the bottom of 
AdS. 

In a typical RHIC experiment, hundreds of nucleons or thousands of quarks are involved 
as shown above. As the interaction is almost simultaneous, thousands of interaction vertices are 
involved and a shower of massive closed strings are created and fall into the center of the AdS space. 
Let N — tt R 2 N Q 2 be the number of such collisions with Q 2 their transverse density and 7r R 2 N the 
transverse nuclcon size. In field theory charged quantas moving with rapidity Y are surrounded 
by extra quanta distributed at smaller rapidities dy = dx/x. In QED the Weiszacker- Williams 
approximation yields a flat distribution of these quantas versus y, i.e. dN/dy constant. In QCD 
dN/dy is not constant and behaves approximately as e a ^^ Y ~ v \ HERA data suggests a s=s 1/4 at 
t w — lGeV 2 . In weak-coupling the BFKL approximation gives ocbfkl{^) = (4a N/ir) In 2, while 
at strong coupling arguments based on AdS/CFT duality yield [16] aAds{t) ~ 7/96 + 0.23 1. In our 
case, we will use a transverse parton density dN/dydx± — Q 2 (y) with Q s in general y-dependent. 

Since the 5th AdS coordinate r is orthogonal to the boundary collision axis, only momenta 
of closed strings with typically Q s are relevant. Then the typical total energy of the closed strings 
is M w J dydN/dyQ s , with all strings assumed to be created instantaneously at the impact. The 
energy of the closed string must be identified as the radial coordinate in the AdS space. The strings 
flake towards the center of the curved AdS space under gravity and arrive at the central region 
simultaneously. The average energy per string is e w Q s . The total energy is therefore of order 
NQ S . We note that N w N 2 and ft « JV C °, so that the total energy is of order N 2 . A black-hole 
forms when the horizon radius is bigger than the size of the closed string. For the large black hole, 
the horizon distance is rt,h = (R^G^M) 1 '' 1 , where G5 is the 5 dimensional Newton's constant. Hence 
the black hole formation condition is 

Q; 1 < (i? 2 G 5 M) 1 / 4 . (46) 

In terms of N(s) the number of pair collisions at the boundary, (46) reads 

N > R- 2 G b Q-\ (47) 

We recall that M w Ql and N w Q 2 rj s a . The entropy generated at the surface by the RHIC 
collision depends on the collision energy as follows 

S^T 3 V 3 ^s 9a ^ 2 (48) 

or an initial temperature T w s 9 "/ 8 . 

D A correction to the big bang picture 

The situation is schematically illustrated in Fig. 7 (b). 

The temperature at the center of the fire ball is given by the warping factor as determined 
above using the distance r from the black hole center to the fire ball center. The temperature at all 
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Figure 7: Sketches of the brane motion for the Big Bang (a) and Little bang (b) geometries, (a) refers 
to a brane moving away from large "black brane" with a time-dependent distance in 5-th dimmension 
r(t). (b) refers to an asymptotically flat brane with the brane to black hole distance r(t), while p(t) is the 
effective spherical size of the fireball on the brane. The Howking radiation has homogeneous time-dependent 
temperature in (a), while it depends on p in case (b). 

other points in the moving brane is warped further since the effective distance now is \/r±(t) 2 + p 2 , 
and vanishes asymptotically. While the precise determination of the imbedding of the moving brane 
is very involved, its shape in the region far from the black hole can be readily discussed approximately 
by neglecting the bending effect. More precisely, from the metric 

ds 2 = -fdt 2 + f^dr 2 + r 2 dn 2 , with / = 1 - G 5 M/r 2 + (r/R) 2 , (49) 

the temperature is approximately given by 

T(t) - Tov/JM/v/l - G 5 M/(r ± (t)2 + + (r ± (t) 2 + p^)/R\ (50) 

where T is the black hole temperature at a reference distance say r = r . Now, the temperature 
depends both on the spatial size p of the fire ball on the boundary, and the distance to the black hole 
r_i_(t). This is the temperature profile which has a peak at the fireball center, p = 0, and decreases 
as p increases. See figure 8. The upshot of this analysis is that the presence of a distance black 
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Figure 8: Plot ofthe temperature for < p < 15 as the brane position r±(t) moves from 5 to 15. The 
bending of the brane is ignored assuming that the brane is far from the black hole. We have used G$M = 1, 
R — 1, Toy/ f(ro) = 1. Notice that the warping factor (~ r 2 ), which will overcompensate the apparent 
shrinking in the fireball size, is not taken into account in the plot. 

hole in bulk produces small additional forces on moving matter in the moving brane. Even without 
considering the bending effect of the brane, the presence of blackhole is detected through the metric. 
For instance, hydrodynamical flow of matter on the brane is now described by T^. v = 0. The source 
of the Christoffels T are two fold: one is the expansion induced by the motion of the brane inside a 
warped background, the other is due to the presence of the black hole. These modifications to ideal 
hydrodynamics are not small even at late stages as far as the strong character of the interaction 
sustains. However, these effects are small far from the fireball center at all times. These analysis 
will provide yet another route to non-ideal hydrodynamics for jV=4 SYM theory at strong coupling. 
In particular to the calculation of the transport coefficients. This will be reported elsewhere. 
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